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We investigate the nature of the interaction-driven Mott-Hubbard transition of the half-filled 
t\—t% Hubbard model in one dimension, using a full-fledged variational Monte Carlo approach 
including a distance-dependent Jastrow factor and backflow correlations. We present data for the 
evolution of the magnetic properties across the Mott-Hubbard transition and on the commensurate 
to incommensurate transition in the insulating state. Analyzing renormalized excitation spectra, we 
find that the Fermi surface renormalizes to perfect nesting right at the Mott-Hubbard transition in 
the insulating state, with a first-order reorganization when crossing into the conducting state. 
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I. INTRODUCTION 



II. MODEL 



Low-dimensional Fermi gases and insulators with in- 
termediate and strong couplings show a plethora of in- 
teresting phenomena, both in the domains of synthesiz- 
able materials 1 and of ultra-cold atom gases, 2 with the 
proximity of metallic, magnetic, superconducting, and in- 
sulating phases being a key target for experimental and 
theoretical studies. One-dimensional correlated electron 
systems are hence good targets, to give an example, for 
the exploration of photo-induced phase transitions^ hav- 
ing in part extremely large third-order non-linear optical 
susceptibilities, with possible applications to all-optical 
switching devices^ 

Here, we are interested in the nature of the interaction- 
driven Mott-Hubbard transition which occurs in the one- 
dimensional t\— 1% Hubbard model at half filling. In par- 
ticular, we will assess the evolution of the Fermi surface 
by varying the Coulomb interaction. In the insulating 
state, the underlying Fermi-surface is given by the bound- 
ary of the occupied states of the renormalized dispersion 
relation, when the residual interactions giving rise to the 
charge gap are turned off in a Gcdanken experiment^— 
Mathematically, the underlying Fermi surface is defined 
in a non Fermi-liquid state as the locus in k— space where 
the real part of the one-particle Green's function changes 
its signi^ii^ By investigating magnetic and charge prop- 
erties, we find that the Fermi surface reconstructs in a 
first-order manner right at the Mott transition. In par- 
ticular, the Fermi surface is generic, namely non-nesting, 
in the metallic side, whereas it has perfect nesting prop- 
erties in the insulating state, at the transition point. 

The paper is organized as follow: in section [Hi we in- 
troduce the Hamiltonian; in section Hill we describe our 
variational wave function; in section IIV1 we present our 
numerical results and, finally, in section [V] we draw the 
conclusions. 



We consider the one-dimensional £i— £2 Hubbard model 
H = - *n4* c i+n,<r + H - c - + ^" i,t " a ' (1) 

i,er,n— 1,2 i 

where c\ a is the electron creation operator, a =f, I the 

electron spin, i — 1, . . . ,L the site index, rii t „ — c\ rT c i a 
the electron density, t\ and £2 the nearest and next- 
nearest neighbor hopping amplitudes,— and U the on-site 
Coulomb repulsion. In this work, we focus our attention 
on the half-filled case with L electrons on L sites. 

The ground state of the t\— £2 Hubbard model at half 
filling is predicted to be an insulator with gapless spin ex- 
citations (conventionally labeled as C0S1) for tijt\ < 1/2 
and every finite U /iir^ a spin-gapped metal (C1S0) with 
strong superconducting fluctuations for £ 2 / > 1/2 and 
small U/tir^ and a fully-gapped spontaneously dimer- 
ized insulator (COSO) for t 2 /h > 1/2 and large U/tiM^ 
Our findings, which are summarized in Fig. [I] are in 
very good agreement with these results. The locus 
of the metal-insulator transition has been investigated 
by several groups ^r— with slightly varying outcomes. 
Remarkably, a transition between incommensurate and 
commensurate spin excitations is expected to take place 
inside the COSO phaseJ^^^ Finally, we would like to 
mention that a tiny C2S2 phase could be stable for 
U/ti — > 0, as suggested by a weak-coupling renormaliza- 
tion group approach^ recent calculations showed that 
this phase can be further stabilized in presence of long- 
range interactions^ 



III. VARIATIONAL APPROACH 

In this paper, we present a variational Monte Carlo 
study of the Hubbard model for £2/^1 > 1/2, which al- 
lows us to determine accurately the locus of the metal- 
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FIG. 1: (Color online) Phase diagram of the t l ~t2 Hubbard 
model at half-filling with the metallic phase with gapped spin 
excitations (C1S0) and the insulating phase with gapless spin 
excitations (C0S1). The insulating phase with gapped spin 
excitations for larger U and £2/^1 > 1/2 has regions with 
commensurate (Q = n) and incommensurate (Q incomm) 
spin-spin correlations. A crossover region separates the phase 
where the peak in S(q) is incommensurate and the one with 
the peak commensurate to a doubled unit cell (Q ~ f/2). 



insulator transition, to study the transition between com- 
mensurate and incommensurate spin-spin correlations 
in the large-C/ (dimerized) phase and to investigate its 
underlying Fermi surface. In particular, we will show 
that the magnetic correlations are related to the single- 
particle spectrum in the optimized variational wave func- 
tion. Moreover, we will propose that the metal-insulator 
transition is driven in the Mott state by a renormaliza- 
tion of the underlying Fermi surface to perfect nesting. 

Both the metallic and the insulating phases can be 
constructed, in a variational approach. In a first step, 
one constructs uncorrelated wave functions given by 
the ground state |BCS) of a superconducting Bardeen- 
Cooper-Schrieffer (BCS) Hamiltonian ; 23 ! 24 

hbcs - + E A i c lt c - q ,i + Hc -> ( 2 ) 

q,a q 

where both the free-band dispersion e q and the pairing 
amplitudes A q are variational functions. We use the 
parametrization 

e q — — 2ti coaq — 2i 2 cos(2q) — fi , 
A q = Aicosq + A 2 cos(2(7) + A 3 cos(3q), ' ' 

where the effective hopping amplitudes t\ and ti 1 as well 
as the effective chemical potential /i and the local pair- 
ing fields Ai, A2, and A3 are variational parameters to 
be optimized. The excitation spectrum for Bogoliubov 
excitations is given by 



U/t 1 =4,t 2 /t 1 =0.75 □ 




FIG. 2: (Color online) For a chain with L — 120 sites, 
the density-density correlations N(q), divided by the mo- 
mentum q, across the metal-insulator transition for i2/ti = 
0.75, 1.1, 1.5. The metallic (insulating) state is characterized 
by a finite (vanishing) value of N(q)/q, in the limit q — > 0. 



The correlated state I^bcs) is then given by 



I^bcs) = J1BCS), 



(5) 



E„ 



(4) 



where J = exp(— 1/2 ■ VijUiUj) is a density-density 
Jastrow factor (including the on-site Gutzwiller term), 
with the Vij being optimized independently for ever dis- 
tance \i — j\. Notably, within this kind of wave function, 
it is possible to obtain a pure (i.e., non- magnetic) Mott 
insulator by considering a sufficiently strong Jastrow fac- 
tor,— i.e., v q ~ l/q 2 (v q being the Fourier transform of 
Vij) and a Luttinger-liquid wave function with arbitrary 
critical exponents^ whenever Vij ~ log \i — j\. In addi- 
tion, a dimerized phase can be obtained just by consider- 
ing a gapped BCS spectrum E q together with v q ~ l/q 2 
(that is the case whenever t 2 /ti > 1/2 and U/ti is large 
enough). 25 Remarkably, in this case, finite dimer-dimer 
correlations are found at large distance, even though the 
wave function does not break the translational symme- 
try. Here, we do not report results on dimer-dimer cor- 
relations, that are found in the C0S0 phase (see Ref . |25[) . 
but we concentrate on spin and charge properties, with a 
particular emphasis on the evolution of the Fermi surface 
by changing ^Ai and U/t\. 

As we demonstrated recentlyj21 the projected BCS 
state I^bcs) can be improved further by considering 
backfiow correlations, which modify the single-particle 
orbitals, in the same spirit as proposed by Feynman and 
Cohcni^S In this way, already the determinant part of 
the wave function includes now correlation effects. All 
results presented here are obtained by fully incorporat- 
ing the backfiow corrections and optimizing individually 
every variational parameter in e q and A,j of Eq. ([3]), in 
the Jastrow factor J of Eq. ([5]), as well as backfiow cor- 
rections. 
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IV. RESULTS 
A. Mott-Hubbard transition 

The ground-state properties can be easily assessed by 
computing density and magnetic structure factors: 

N(q) = j-^e^ k - l Hn k ni), (6) 
fc,; 

S(q) = ~$>^)<S|Sf>, (7) 
fc,; 

where and Sf, are the total density and the z- 
component of the spin operator on site k, respectively. 

The static density-density correlations behave qualita- 
tively different in a metallic and a Mott-insulating state 
for small momenta q, with the metallic state being char- 
acterized by a linear dependence of N(q) ~ q, while in 
the insulating phase N(q) ~ q 2 2^ In Fig. [5J we present 
the behavior of N(q)/q across the transition for three 
values of the ratio t 2 jt\. The locus of the Mott-Hubbard 
transition can be determined easily, allowing us to draw 
the phase diagram in Fig. [T] Our determination of the 
line separating the metallic and the insulating phase is 
in good agreement with Refs. [ililfLBl 

The metallic region in the phase diagram can be de- 
scribed as a Luther-Emery liquid, with a finite gap in 
the spin excitation spectrum and gapless charge excita- 
tions^ The charge stiffness K p can be extracted, for 
example, from the long-distance behavior of the density- 
density correlations. In any conducting phase, we expect 
that K p is also related to the slope of N{q) at small q, i.e., 
N(q) ~ K p \q\/ir. In fact, the latter equation, which is 
definitely valid in Luttinger liquids, should hold whenever 
the charge degrees of freedom are gapless^ This proce- 
dure to obtain K p works very well for the doped single- 
band Hubbard model, namely for the model of Eq. ([1]) 
with t 2 = 0f^2 in comparison with the exact results, as 
obtained by Bethe ansatz.— Unfortunately, in the metal- 
lic phase of the ti—t 2 model, the existence of very small 
gaps prevents us to obtain accurate results for K p . 

B. Magnetic properties 

The presence of short-range magnetic order is signaled 
by the appearance of a peak in S(q), for a certain momen- 
tum Q. In the following, we will compare the magnetic 
properties with the renormalized single-particle spectrum 
E q of the optimized variational wave function. The en- 
ergy scale for E q will be taken as the bandwidth W of 
the original free dispersion e q = — 2ti cos(q) — 2t 2 cos(2<7). 
Note that, in the non- interacting case, there is only a sin- 
gle, perfectly nested Fermi surface for t 2 /ti < 0.5, with 
two Fermi points separated by n. Instead, for t 2 /ti > 0.5 
there are two Fermi seas and four Fermi points. 

In the metallic phase, the spin properties are only 
slightly modified by the presence of a small but finite 




FIG. 3: (Color online) Upper panel: Spin-spin correlations 
S(q) at £2/^1 = 0.75, for a L — 120 chain. Data are shown for 
U/ti =0,4 (metal) and for U/ti — 5,6, 12 (insulator). Lower 
panel: Single-particle spectrum E q /W at t^jt\ — 0.75 for the 
same values of the electron-electron repulsion U and the same 
chain length. 



interaction U, with respect to the U = behavior; for 
t 2 /ti > 0.5 the single-particle spectrum E q exhibits four 
minima, at ±fci and ik 2 , and the peak of S(q) is located 
at Q mct = k 2 - fci = tt/2. The condition Q mot = tt/2 
is determined by the Luttinger sum-rule for the metal, 
which states that the total volume of the Fermi sea equals 
the number of electrons. In the insulating phase, the sit- 
uation changes qualitatively and the magnetic properties 
of the system become strongly affected by the value of 
ta/h. 

In Fig. [3J we show the behavior of S(q) across the 
metal-insulator transition for t 2 /ti = 0.75, in comparison 
with the variationally determined renormalized single- 
particle spectrum E q . It can be observed that, when 
entering the insulating phase, the single-particle spec- 
trum becomes strongly gapped and the two central min- 
ima collapse into a unique relative minimum at q = 0, 
that subsequently disappears, as U /t\ increases. At the 
same time, the peak in S(q) shifts from Q mct = tt/2 
to Q ms = it. Remarkably, just above the Mott transi- 
tion, namely for 5 < U/ti < 10, the quantity 2k\ (i.e., 
the distance between the two absolute minima of E q ) 
is slightly different from 7r and becomes commensurate 
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FIG. 4: (Color online) Evolution of 2fci (distance between 
the minima of the single-particle spectrum E q ) in the insu- 
lating phase for £2/^1 = 0.75. Note the transition from an 
incommensurate to a commensurate value for U/ti ~ 9. 



only after a second transition (e.g., U/ti ~ 9), inside the 
insulating phased see Fig. |U However, the degree of 
incommensurability is very small and does not show up 
in corresponding shift of the peak in S(q) from Q ms = n. 
Indeed, the magnetic correlations are short-ranged and 
the peak in S(q) is consequently broad. A shift of the 
momentum away from tt will therefore result in a shift 
of the maximum in S(q) only for a substantial degree of 
incommensurability. 

In Fig. 03 we plot the spin-spin correlations S(q) and 
the single-particle spectrum E q for the ratio ia/ii = 0.9. 
In the metallic phase, the spin-spin correlations are al- 
ways peaked at Q mct = ir/2, while in the insulating phase 
the peak slowly shifts to Q — O.671-. For U = 6 and 7 the 
single-particle spectrum E q is qualitatively different from 
the one for larger J7's. As shown later, this is related to 
the different behavior of the variational hopping ratio 
til~k\ close to the metal-insulator transition with respect 
to the strong-coupling regime. For larger values of the 
ratio U/t±, E q shows four local minima, with the peak 
in S{q) located at Q lns — 2fci, where 2k\ is the distance 
between the two absolute minima. 

Finally, in Fig. [5J we summarize the spin-spin correla- 
tions for different t%/t\, at a given value of U/ti, chosen 
to be far enough from the metal-insulator transition in 
order to describe the large-C/ behavior of S(q). The peak 
in the spin-spin correlations exhibits the commensurate- 
incommensurate transition moving far from Q = n, as 
the ratio ti/t\ is increased. When ti/t\ = 1.5 the system 
behaves already like in the £2/^1 ~^ 00 limit, with the 
peak commensurate to a lattice with a doubled unit cell 
(Q = tt/2). These results are in agreement with previous 
studies for the Heisenberg^ and the Hubbard models 
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FIG. 5: (Color online) Upper panel: Spin-spin correlations 
S(q) at £2/^1 = 0.9, for a L = 120 chain. Data are shown 
for U/ti — 5 (metal) and U/ti = 6, 7, 8, 12 (insulator). Lower 
panel: Single-particle spectrum E q /W at t%/ti = 0.9 for the 
same values of the electron-electron repulsion U and the same 
chain length. 



C. Fermi surface renormalization 

Finally, we present our central result, namely the 
fact that the metal-insulator transition is driven, in the 
Mott-insulating state, by a renormalization of the un- 
derlying Fermi surface to perfect nesting. With under- 
lying Fermi surface, we mean the locus of the high- 
est occupied momenta in the non-interacting spectrum 
e g = —2tiCOs(q) — 2i2Cos(2g), obtained from the opti- 
mized variational hopping parameters. We would like to 
stress that the concept of an underlying Fermi surface is 
of central importance for the angular resolved photoemis- 
sion spectroscopy (ARPES) studies of strongly correlated 
systems, like the high-temperature superconductors^— 

Note, that E q = e q + corresponds within renor- 

malized mean- field theory 2 ^ to the excitation spectrum of 
projected Bogoliubov quasiparticles and e q hence to the 
dispersion of the renormalized quasiparticles. Moreover, 
recent calculations on the t—J and the periodic Ander- 
son models highlighted the possibility to assess the Fermi 
surface from the parametrization of a variational wave 
function^i^i Here, the renormalization of the hopping 
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FIG. 6: (Color online) The spin-spin correlations S(q) at 
U/ti = 12 and for L = 120 sites, for different values of the 
hopping ratio tilt\. Arrows indicate the quantity 2fci, ob- 
tained from the single-particle spectrum E q . 

parameters made it possible to show non-trivial defor- 
mations of the non-interacting Fermi surface, due to the 
Gutzwiller projection. 

We show in Fig. [7] that the ratio f 2 /*i m the metallic 
phase is almost equal to the bare value t%jt\, regard- 
less of the degree of interaction. This weak renormal- 
ization of the band-structure in the metallic state is in 
agreement with a renormalization-group study, 21 which 
predicts that the renormalization of the Fermi surface is 
proportional to U 2 . Then, after the metal- insulator tran- 
sition, the ratio jumps to a smaller value, very close to 
1/2. According to our data, we propose that the opti- 
mized variational ratio of t%/ti is renormalized to 1/2, 
exactly at the metal-insulator transition. This discontin- 
uous behavior of the renormalized band structure is also 
evident in Fig.[5j e.g., for t 2 /ti = 0.9, with the number of 
minima of the single-particle spectrum E q jumping from 
four to two when entering the Mott-insulating state. We 
note that an analogous tendency towards a Fermi surface 
symmetrization in the insulating state has been observed 
in a study of a two-dimensional frustrated lattice! 35 ' 36 

A renormalization of the variational hopping ratio to 
ta/ti = 1/2 implies that the Fermi surface is nested, with 
two Fermi points separated by a vector ir. This perfect 
nesting condition drives the system to be an insulator, 
generating a charge gap as soon as electron-electron in- 
teraction is switched on. Remarkably, while the metal- 
insulator transition is driven by the renormalized disper- 
sion e q , the pairing terms A q are crucial in determining 
the spin properties of the model, via the renormalized 

excitation spectra E q — J e 2 + A 2 . Indeed, as shown for 

example in Fig.[5j the minima of the single-particle spec- 
trum at t2/t\ = 0.9 are connected by an incommensurate 



vector, leading to an incommensurate peak in S(q). 
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FIG. 7: (Color online) The variationally optimized hopping 
ratio t2/t\ in |BCS), see Eq. as a function of U/t\. The 
metal-insulator transition takes place for U/ti = 4.5 ± 0.5, 
5.5±0.5 and 6.5±0.5 for t 2 /ii = 0.75 (triangles), 0.9 (squares) 
and 1.0 (circles), respectively. The variational state jBCS) 
contains only a single Fermi sea for t2/t\ < 0.5 (horizontal 
line). 

V. CONCLUSIONS 

We have presented an extensive study of the phase 
diagram of the one-dimensional t\— 1% Hubbard model 
at half filling, with emphasis on the evolution of the 
magnetic properties and of the underlying Fermi surface 
across the interaction-driven Mott-Hubbard transition. 
We have shown that the magnetic correlations are re- 
lated to the single-particle spectrum in the optimized 
variational wave function and we have described how 
they are affected by the metal-insulator transition. In 
the insulating phase, the peak in the spin-spin correla- 
tions exhibits the commensurate-incommensurate tran- 
sition moving far from Q = 7r, as the ratio ^Ai is in- 
creased, and then becomes commensurate to a doubled 
unit cell when ^Ai ^ 1-3. 

Our main findings culminate in the hypothesis that the 
underlying Fermi surface renormalizes to perfect nesting 
right at the transition in the insulating phase, with a first- 
order reorganization when crossing the transition into the 
metallic state. Similar renormalizations of the Fermi sur- 
face have been observed in two-dimensional models^ ' 35 ' 36 
Therefore, we believe that our results are important 
for an improved understanding of Mott-Hubbard tran- 
sitions quite in general, transcending the specific one- 
dimensional physics. 

L.F.T. and C.G. acknowledge the support of the Ger- 
man Science Foundation through the Transregio 49. 
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